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— - We show that a system of polaritons - combined atom and photon excitations - in an 

array of coupled cavities, under an experimental set-up usually considered in electromag- 
netically induced transparency, is described by the Bose-Hubbard model. This opens up 
the possibility of using this system as a quantum simulator, allowing for the observation 
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of quantum phase transitions and for the measurement of local properties, such as single 
site observables. All the basic building blocks of the proposed setting have already been 
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achieved experimentally, showing the feasibility of its realization in the near future. 
1. INTRODUCTION 



> 

The exponential growth of the Hilbert space dimension with the number of particles of the 
system is a hurdle for the classical simulation of quantum many-body systems, in particular near 
or at criticality. A possible solution to this conundrum, pioneered by Feynman 20 years ago HI, 
is the use of quantum simulators, i.e. artificial quantum systems which generate the evolution of 
rather complex Hamiltonians in a controlled and systematic manner. An emblematic example 
of a quantum simulator is a quantum computer, capable of simulating any quantum evolution. 
Yet, the overwhelmingly high degree of accuracy required to build a useful quantum computer - 
operating below the fault-tolerance threshold - renders unlikely the availability of such a device 
in the near future. A natural approach is then to consider less general quantum simulators, engi- 
neered to simulate specific many-body models, but requiring much less stringent experimental 
conditions. 

Cold atoms trapped in an optical lattice constitutes an important example of this second 



kind of simulator [|2l[3l|4l[5]|. By adjusting the experimental parameters properly, the Bose- 
Hubbard (BH) Hamiltonian - extensively used in modeling strongly correlated many-particle 
systems - can be generated. This idea, in turn, culminated in several important experiments em- 
ulating condensed-matter phenomena, such as the superfiuid-to-Mott-insulator quantum phase 
transition [3] and the generation of a Tonks-Girardeau gas for the first time 0. Moreover, the 
realization of the BH model opened up the possibility of studying several interesting aspects 
of many-body quantum systems, such as disorder, Anderson localization, and quantum mag- 
netism ||5). Despite their impressive success, optical lattices are in all their applications limited 
by an intrinsic draw back: due to the separation between neighboring lattice sites of only half 
of the employed optical wavelength, it is extremely challenging to access them individually. 

Recently, we have proposed an alternative way to create an effective Bose-Hubbard model, 
which, in contrast to optical lattices, allows for the manipulation and measurement of the prop- 
erties of individual constituent particles H (see Refs. Q for related proposals). This proposal 
can be realized in coupled arrays of optical cavities each interacting with an ensemble of atoms 



Figure 1: (taken from Ref. [6]) A schematic view of the model considered. 



leading to the dynamical evolution of polaritons, combined atom-photon excitations, in this 
setting. An exciting feature of this proposal is the possibility of observing strong-correlation 
phenomena in a system of photons, which emerges as one of the limits of the polariton excita- 
tions. 

Our model consists of an array of cavities in an arbitrary geometry, where each cavity inter- 
acts with an ensemble of atoms, which are driven by an external laser. Photon hopping occurs 
between neighboring cavities due to the overlapp of their photonic wavefunctions (see fig. [Q) 
while the repulsive force between two polaritons occupying the same cavity is generated by a 
large Kerr nonlinearity that occurs if atoms with a specific level structure, usually considered in 
the context of Electromagnetically Induced Transparency [8 , 9 ll0lllllll2| , interact with light. By 
varying the intensity of the driving laser and hence the strength of the generated Kerr nonlinear- 
ity, the system can be driven through the superfluid-to-Mott-insulator transition. In particular, 
the driving laser may be adjusted for each cavity individually allowing for a much wider range 
of tuning possibilities than in an optical lattice, including attractive interactions 10. 

Our system can be described in terms of three species of polaritons, where under the con- 
ditions specified below, the species which is least vulnerable to decay processes is governed by 
the effective Bose-Hubbard (BH) Hamiltonian, 

#eff = * E (4) 2 w 2 + J E (pW + h - c -) • (1) 

R (R,R') 

creates a polariton in the cavity at site R and the parameters k and J describe on site repul- 

R 

sion and inter cavity hopping respectively. 

2. DERIVATION OF THE MODEL 

We start by considering a periodic array of cavities, which we describe here by a real, 
periodic dielectric constant. The electromagnetic field can be expanded in Wannier functions, 
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Figure 2: (taken from Ref. 0) Diagram of the level structure considered. 

vug, each localized at one single cavity at location R. In terms of the creation and annihilation 
operators of the Wannier modes, at and as, the Hamiltonian of the field can be written, 

Ft 

n = u c ^2 \ a n a R + \ ) + 2uJca ( a W' + h - c -) • <2) 

R (R,R r ) 

Here is me sum °f ai l P arrs °f cavities which are nearest neighbors of each other, a is 

given by an overlap integral of neighboring Wannier functions [fT3l and can be obtained numer- 
ically for specific models. Since a < 1, we neglected rotating terms which contain products 
of two creation or two annihilation operators of Wannier modes in deriving ©. The model 
© provides an excellent approximation to many relevant implementations such as coupled 
photonic crystal micro-cavities or fiber coupled toroidal micro-cavities. Furthermore, it allows 
for the observation of state transfer and entanglement dynamics and propagation for Gaussian 
states lfl4l[Bl 

The on site interaction potential for polaritons is generated in each cavity by atoms of a 
particular level structure that are driven with an external laser in the same manner as in Elec- 
tromagnetically Induced Transparency, see figure 2: the transitions between levels 2-4 and 1-3 
couple via dipole moments to the cavity resonance mode, whereas the transitions between levels 
2 and 3 are coupled to the laser field. 

In a rotating frame with respect to H = ujc {a) a + |) + J2f=i (^0^22 + ^033 + ^c^m) > 
the Hamiltonian of the atoms in one cavity reads, Hj = ^=1 ( £cr 22 + ^°33 + (A + e)a{^) + 
Yl!j=i (^i a 23 + #i3°i3 a ^ + 024 024 a* + h.c.) . Here a\, x = \kj)(lj\ projects level / of atom 
j to level k of the same atom, a) creates one photon in the cavity, ujc is the frequency of the 
cavity mode, Vt l the Rabi frequency of the driving by the laser, gi 3 and g 2 4 the parameters of 
the dipole coupling of the cavity mode to the respective atomic transitions and A, S and e are 
detunings. All these parameters are assumed to be real. 



Neglecting the coupling to level 4, g 2 4 = 0, and two photon detuning, e — 0, the Hamil- 
tonian Hj can be written as a model of three species of polaritons 11 1 Oil 1 611 with the creation (and 

annihilation) operators pi = ^ (gS\ 2 — fi^a 1 ^ and p± = A ^±s) {^lS\ 2 + W ± ^S\ 3 

where g = y/N 9l3 , B = y^+W L , A = V45 2 + 6 2 , S\ 2 = ^Ef=i4 and S\ 3 = 

"7KF ^i=i °3i' ^ n tne °^ l ar S e atom numbers, N ^> 1, pj,, p+ and pi satisfy bosonic com- 
mutation relations and for g 2 4 = and e = the Hamiltonian if/ thus describes independent 
non-interacting bosonic particles, 

[#/] S24=0i£=0 = no pIpo + fi+ p + p+ + n- pip- , (3) 

where the polariton frequencies are given by p,Q = 0, /i + = (6 — A)/2 and /i_ = (5 + A)/2. 
We will now focus on the dark state polaritons, p\, which will be shown to be described by the 
Hamiltonian (HJ). 

For \g 2 A , \e\, |A| <C — /Uq| , — /x |, the coupling g 2 4 and the two photon detuning 
£ do not induce interactions between different polariton species [6J. Assuming further that 
\Jn p (n p — 1)^24 I A|, where n p is the number of polaritons p\, the coupling to level 4 can be 
treated perturbatively. The result is an energy shift of n p {n p — 1)k with 

K= _gk m^i 2 . (4) 

a (jv 9 ? 3 + ny 2 

Note that k > for A < and vice versa. The two photon detuning e, in turn, leads analogously 
to an energy shift of e g 2 B~ 2 for the polariton p\, which plays the role of a chemical potential 
in the effective Hamiltonian. 

The interaction between cavities, derived above in terms of the photonic operators, can be 
reexpressed in terms of the polaritons scpecies. Due to the large separation in the polaritons 
frequencies, |2u;<7a| <C — fi \ , \fi- — fi \, the hopping does not mix different polariton 
species. The Hamiltonian for the polaritons (the dark state polariton at site R) thus takes 
on the form (OQ), with 

where we have assumed a negligible two photon detuning, e ~ 0. 

The number of polaritons in one individual cavity can be measured via resonance fluores- 
cence. To that end, the polaritons are made purely atomic excitations by adiabatically switching 
off the driving laser O. The process is adiabatic if gB~ 2 j- t Vl L <C |//_| H21, which 
means it can be done fast enough to prevent polaritons from hopping between cavities during 
the switching. Hence, in each cavity, the final number of atomic excitations in level 2 is equal 
to the initial number of dark state polaritons. By measuring the number of polaritons in several 
individual cavities in possibly repeated resonance flourescence experiments one can obtain the 
local number fluctuations which allow to distinguish between the different phases of the system. 



3. FEASIBILITY OF THE OBSERVATION OF PHASE TRANSITIONS 



Promising candidates for an experimental realisation are photonic band gap cavities H171I18H 
and toroidal or spherical micro-cavities, which are coupled via tapered optical fibres [[P9ll . These 



cavities can be produced and positioned with high precision and in large numbers. They have 
a very large Q-factor (> 10 8 ) for light that is trapped as whispering gallery modes and efficient 
coupling to optical fibres [20] as well as coupling to Cs-atoms in close proximity to the cavity 
via the evanescent field II2TII221 have been demonstrated experimentally. Photonic crystals 
represent an appealing alternative as they offer the possibility for the fabrication of large arrays 
of cavities in lattices or networks lfT3ll23ll . These cavities have been realised with Q-factors of 
10 6 and higher Q-factors of 2 x 10 7 have been predicted ||T8l . 

As levels 1 and 2 are metastable, spontaneous emission can only occur from levels 3 and 4. 
However, due to pratically negligible population of these two levels, the main loss channel is 
cavity decay. For successfully observing the dynamics and phases of the effective Hamiltonian 
(1), the repulsion term k needs to be much larger than the damping rate for the dark state 
polaritons, T [6]. Choosing CIl <C one can currently achieve a ratio of k/T ~ 5.2 for 

photonic band gap cavities (for predicted Q-factors this can go up to k/T ~ 170) IfTBl , while 
for toroidal micro-cavities the impressive large value of k/T ~ 1.1 x 10 3 is possible 112411211 . 
making both ideal candidates for an experimental implementation. 

As an ilustration, a numerical simulation of the Mott-insulator-to-superfluid phase transition 
is shown in figure 3. We consider three coupled cavities with periodic boundary conditions, 
which are prepared such that there is initially one polariton pL in each. We take parameters 
for toroidal microcavities from E4l . g 2 4 = #13 = 2.5 x 10 9 s _1 , T 3 = T 4 = 1.6 x 10 7 s _1 
(spontaneous emittion rate from levels 3 and 4), Tc = 0.4 x 10 5 s _1 (cavity decay rate), N = 
1000, A = -2.0 x 10 10 s~\ and 2u c a = 1.1 x 10 7 s -1 . To drive the system through the 
transition, The Rabi frequency of the driving laser is imcreased from £l L = 7.8 x 10 10 s _1 to 
VIl = 1-1 x 10 12 s _1 throughout the experiment, while all other parameters remain constant. 
Fig. 3 shows in the left part, VLl (inset), the on-site interaction strength k and the hopping rate 
J during the transition. The right part shows that number fluctuations F\ = ((p\pi) 2 ) — (p\pi) 2 
are small in the Mott phase and increase when driving into the superfluid phase. The average 
number of polaritons in one cavity n% = (p\pi) stays close to unity across the whole simulated 
time range confirming that polariton loss is indeed negligible on this time scale. 

4. CONCLUSION 

In summary, we have demonstrated, that polaritons, combined atom - photon excitations, 
can form an effective quantum many particle system, which is described by a Bose-Hubbard 
Hamiltonian. We have demonstrated numerically that the observation of the Mott insulator 
to superfluid phase transition is indeed feasible for experimentally realizable parameters. In 
contrast to earlier realizations, our scenario has the advantage that single sites of the BH model 
can be addressed individually and offers the possibility of creating a BH type Hamiltonian with 
an attractive on site potential. 
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